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We investigate even-even two-proton borromean systems at prominent intermediate heavy waiting 
points for the rapid proton capture process. The most likely single-particle levels are used to calculate 
three-body energy and structure as a function of proton-core resonance energy. We establish a 
linear dependence between two- and three-body energies with the same slope, but the absolute 
value slightly dependent on partial wave structure. Using these relations we estimate low-lying 
excited states in the isotones following the critical waiting points. The capture rate for producing a 
borromean bound state is described based on a full three-body calculation for temperatures about 
0.1 — 10 GK. In addition, a simple rate expression, depending only on a single resonance state, is 
found to comply with the three-body calculation for temperatures between 0.1 and 4 GK. The rate 
calculations are valid for both direct and sequential capture paths. As a result the relevant path 
of the radiative capture reactions can be determined. We present results for El and E 2 photon 
emission, and discuss occurrence preferences in general as well as relative sizes of these most likely 
processes. Finally, we present a method for estimating proton capture rates in the region around 
the critical waiting points. 

PACS numbers: 25.40.Lw, 21.45.-v, 26. 


I. INTRODUCTION 

A number of proton dripline nuclei are of particular in¬ 
terest as waiting points in the rapid-proton capture (rp) 
process expected to be active in the accretion of a close 
binary system containing a neutron star and resulting in 
an X-ray burst W- When the proton binding energy be¬ 
comes negative at the dripline another proton is needed 
to produce the borromean system (a bound system with 
unbound subsystems). The effective lifetime of the critical 
waiting points in a stellar environment is a central quantity 
in the understanding of this astrophysical process. This de¬ 
pends crucially on both proton binding of the waiting point 
plus one and two protons, as well as the reaction rate form¬ 
ing these nuclei 0 ( 3 - The current estimates still result 
in an uncertainty in the effective lifetime of several orders 
of magnitude 00- The energy, capture time and capture 
mechanism can be explored through three-body calcula¬ 
tions. Relatively few full three-body results have been pub¬ 
lished, although many capture rates have been estimated 
using various approximations Enama. So far the three- 
body results have been limited to nuclei lighter than 38 Ca 
mm- However, three heavier critical waiting points ex¬ 
ists, which have as of yet not been treated from a few-body 
perspective. These critical waiting points are 64 Ge, 68 Se, 
and 72 Kr El HI]- Recent efforts suggest 64 Ge is of less im¬ 
portance than previously thought, while 68 Se, on the other 
hand, is thought to be of prime importance 12 13 . 

Weakly bound nuclear states have been successfully de¬ 
scribed as few-body structures in a number of cases for 
a long time [14, [T5J. The most thorough and abundant 
theoretical investigations exist along the neutron dripline 
nn m 3 , but also excited states of ordinary nuclei have re¬ 
vealed this structure 0, and recently a few medium-heavy 


alpha-dripline nuclei were suggested to be of two-alpha + 
core structure [15] . Few-body formalism is most often ap¬ 
plied to the very light nuclei, where the constituent par¬ 
ticles have a less intricate structure. For instance, few- 
body models have been applied successfully in describing 
the bridging of the mass gaps at A = 5 and 8 PHEU]. Like¬ 
wise, three-body models have been applied for two-neutron 
plus core systems to provide astrophysically relevant pro¬ 
duction rates at the neutron dripline mm again for the 
very light systems. The next step to the proton dripline 
nuclei and two-proton structure has been investigated for 
special, relatively light nuclei EJEHEHH], but much less 
for heavier systems. For the heavier nuclei the increased 
Coulomb interaction severely complicates at least the nu¬ 
merical calculations, if not the conceptual picture. 

The location of the proton dripline - which is rather ir¬ 
regular - is fairly well established at least up to around 
the medium-heavy nuclei 0 I24ti26] , but details of the nu¬ 
clear properties are often very scarce. The most important 
quantity is of course the binding energy and the related sta¬ 
bility which provides the definition of the dripline. Near the 
neutron dripline it is established that low-angular momen¬ 
tum single-particle states produce spatially extended two- 
neutron halo nuclei 0- This few-body picture is less effec¬ 
tive at the proton dripline since the unavoidable Coulomb 
interaction would produce a confining barrier for bound 
nuclei. However, a cluster structure of decoupled core and 
proton degrees of freedom still constitutes a fair descrip¬ 
tion provided the core is relatively tightly bound, and the 
few-body binding energy is very small or perhaps even neg¬ 
ative. 

Performing two and three-body calculations at the pro¬ 
ton dripline requires information about the single-proton 
orbits, in particular their energy, as well as spin and parity 


relative to the core nucleus. Measured values are in most 
cases not available for these ground and low-lying excited 
states, which should be occupied by the protons just before 
or after reaching the dripline. Instead, trends from neigh¬ 
bouring nuclei and especially the results from mean-field 
calculations are used. The inherent uncertainties can be 
ascertained by variation of these input parameters. 

The purpose of this paper is first of all to provide charac¬ 
terizing information about the rp process at waiting points 
for intermediate heavy nuclei along the proton dripline. 
The intend is to provide an initial method for approaching 
nuclei around the critical waiting points from a three-body 
perspective. The borromean three-body structure is crucial 
as an intermediate structure which therefore first must be 
investigated. We shall determine structure and constraints 
for the decisive two-body proton-core energies. This is used 
to present a general method for estimating proton capture 
rates given very sparse experimental information. 

In Sec. [IT] we briefly give the theoretical framework and 
pertinent formulae for both the few-body formalism and the 
three-body reaction rates, as well as the parametric choices 
of the potentials. Section |III| contains a detailed descrip¬ 
tion of the ground state properties of the borromean proton 
dripline systems of interest. Both the relation between two- 
and three-body energies and the structure of the ground 
state is examined in detail. Section |IV| gives calculated re¬ 
sults for the radiative capture reaction rates, based on the 
most likely two-body energies. This also includes a neces¬ 
sary examination of the most likely excited states. Section 
0 discusses limiting values of the basic physical parame¬ 
ters, and predict borromean structure and the related reac¬ 
tion mechanism at the astrophysical waiting points. This is 
combined to provide a method for estimating proton cap¬ 
ture rates around the critical waiting points. Finally, we 
briefly summarize and conclude in Sec. m 


II. THEORETICAL FRAMEWORK 


Our aim is to calculate proton dripline three-body struc¬ 
tures and capture rates from continuum two-proton states 
into the corresponding borromean ground state at medium 
heavy waiting points. For this we need properties of ground 
and continuum states as well as the derived rates and cross 
sections. The general theoretical background can be found 
in Refs. mmmm, but for completeness we shall here 
collect the ingredients necessary to explain notation, cal¬ 
culations and results. We shall only discuss two-proton 
three-body systems. 


A. Three body procedure 


We outline briefly our method of hyperspheric adiabatic 
expansion of the Faddeev equations in coordinate space 
EH- The system of two protons and a core can be described 
by three sets of the relative Jacobi coordinates, (x i; yj), de¬ 


fined as 
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where to is a normalization mass chosen to be the nucleon 
mass of 939 MeV/c 2 . The six hyperspherical relative coor¬ 
dinates are the two pairs of directional angles, (fl^, H,^), 
for Xi and y^, and hyperangle, a*, and hyperradial, p, co¬ 
ordinates defined by 


Xi = psin ai, yi = p cos ai . (4) 


The Hamiltonian, H , can be written both in Jacobi and 
hyperspheric relative coordinates. For later convenience we 
focus here on the Jacobi set where x cp connects core and 
proton, that is we can define 
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where p cp and p cPiP are the reduced masses of core-proton 
and proton to core-proton systems, respectively. The three 
two-body interactions are the proton-proton interaction, 
V pp , the core-proton interaction, V cp , and the interaction 
between the core and the second proton, V cp >. The coordi¬ 
nate dependencies of the two-body interactions, V cp (x cp ), 
V cp '(x cp ,y c p), and V pp (x cp , y cp ), between protons and core 
are given as arguments. In hyperspheric coordinates we 
have 
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where l X i and l y i are the angular momentum operators re¬ 
lated to Xj and y 

The method consists of an adiabatic expansion of the 
total wave function 4/, that is 


^ = p- 5/2 J2fn(P)^n(p^) 1 ( 11 ) 

n 

where each of the angular wave functions, is a sum of 
Faddeev components, = </> n ,i+</ ) n, 2 +</ > n, 3 , related to the 
three corresponding Jacobi sets, and obeying the Faddeev 
equations 

Oyri 

0 = (A 2 - A„) + ^p 2 Vi§ n . (12) 
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The angular eigenvalues, A„(p), and the related complete 
set of eigenfunctions, are first computed for each p. 
Subsequently the radial expansion coefficients, f n {p ), are 
found from the coupled set of radial equations 
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The left hand side reveals the crucial effective adiabatic 
part of the diagonal potential acting on the particles 


V, 
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B. Potentials and properties 

The result of the three-body calculations is dictated by 
the two-body potentials employed. The proton-proton in¬ 
teraction in free space is well known in many details and 
with high accuracy. However, its influence on the three- 
body solutions is only marginal provided, first of all, that 
the s-wave scattering length is reproduced, and second that 
the low-energy properties of the p and d partial waves are 
of reasonable (small) size and in fair agreement with the 
experimental values. The phenomenologically adjusted po¬ 
tentials described in Ref. [28] were used for the proton- 
proton interaction. 

The test case used throughout this paper is the three- 
body system 68 Se + p + p ( 70 Kr). However, the results ap¬ 
ply to the region in general, as the small mass and charge 
variations between the three critical waiting points are in¬ 
consequential for our purposes. 

The proton-core potential is on the other hand decisive 
and a careful choice has to be made. For light-to-medium 
heavy cores of mass numbers around 68 we use the Woods- 
Saxon form with a spin-orbit interaction, that is 


no=r c(r ) + l+e( ;_ w . 

+ i.slA , Xv • m 

r dr 1 + eb~ R is)! a is 

where r is the distance between the two bodies, Vc is the 
Coulomb potential, V 0 and Vq s are the potential strengths 
of the nuclear and spin-orbit potential, 1 and s are the angu¬ 
lar momentum and spin operators, and R and a are param¬ 
eters governing the radius and the thickness of the poten¬ 
tials. This form is used for all partial waves, although the 
two strength parameters vary strongly depending mainly 
on the angular momentum. 


A three-body potential will generally not be included. 
This would eventually be needed to adjust the energy lev¬ 
els according to a measured two- and three-body energy 
spectrum. However, the effect of such an addition is con¬ 
sidered throughout the paper. 

For the mass region in question we choose the values of 
the radial shape parameter to be R = 7.2 fm, Ri s = 6.3 fm, 
a = 0.65 fm, and ai s = 0.5 fm. These choices are motivated 
by the knowledge of the average nuclear mean-field poten¬ 
tials and densities [33]. Accurate values are not needed 
because adjustments of the strengths in any case are neces¬ 
sary for fine-tuning of the two-body energies. This is also 
one reason for omitting more complicated spin-dependence 
like the tensor or quadratic spin-orbit potentials. In addi¬ 
tion, the Coulomb potential is chosen from homogeneous 
charge distributions of radii, 5.6 fm, and 1.8 fm for core¬ 
proton system and proton-proton systems, respectively. 

The two-body strength parameters are usually adjusted 
to reproduce the core-proton bound and resonance energies, 
but unfortunately the energy spectrum of 69 Br is not known 
experimentally. Based on shell model calculations the re¬ 
gion is known to be near the midpoint of the fpg- shell. The 
most likely two-body orbitals are then f 5 / 2 , p 3 / 2 (or possi¬ 
bly P 1 / 2 ), and gy/ 2 - However, recent experiment show that 
the gy / 2 orbital is not important for nuclei around A = 70 
with TV ~ Z (29]. Instead / and p orbitals are assumed to 
dominate the low-lying spectrum. This is also confirmed 
by the known spectrum of the mirror nuclei [i3U • To get 
opposite parity single-particle orbitals, necessary to form 
negative parity three-body states, a d $/2 orbital can be in¬ 
cluded instead. To allow occupancy only of these selected 
two-body core-proton states in the three-body calculation is 
tricky because both lower- and higher-lying levels must be 
excluded. The Pauli forbidden two-body bound states are 
excluded for each partial wave by use of shallow potentials 
without bound states. The large-distance properties are 
then precisely correct but the unimportant nodes at small 
distances are then not present for these excited states. 

To locate one and only one partial wave at a given small 
energy we must provide an accurately adjusted attractive 
potential, while all other partial waves must have suffi¬ 
ciently strongly repulsive potentials to prevent occupancy. 
A given two-body energy then provides one constraint cor¬ 
relating the two strengths. To select one and only one of 
two spin-orbit partners we choose a relatively high, positive 
or negative, value of Vy s . This may result in an abnormal 
order of spin-orbit partners, but the goal is achieved. The 
different partial waves are completely independent of each 
other on the two-body level, and we can therefore place all 
of them as we choose. 


C. Radiative capture rate 

We want to calculate the waiting time before two pro¬ 
tons from the astrophysical environment are captured by 
a proton-dripline nuclear bound core. The reaction rate, 
i?, for the corresponding one-step 7, three-body transition 
process, p + p + c -t A + 7, is given in general in Ref. [3Tj. 
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For the special system of two protons plus an even-even 
core the rate becomes 
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where E 1 = E + B is the photon energy, E the total three- 
body energy, and B is the three-body (positive) binding 
energy of the even-even nucleus, A, with the wave function, 
'F 0 . The photodissociation cross section, <t 7 (.E 7 ), for the 
inverse process A +7 —>■ a + b+c, is a sum over contributing 
electric multipole transitions of different orders, £. That is 
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where the strength function for the ££ transition, 


dE 
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is given 
(,(*) 


through the reduced matrix elements, 
(^ tJ || 0 <||\f r o)j where 0 ^ is the electric multipole op¬ 
erator, ij/f 1 is the wave function of energy, Ei , for all 
bound and (discretized) three-body continuum states in 
the summation. 

The astrophysical processes most often occur in a gas of 
given temperature, T, which means we have to average the 
rate in Eq. (18) over the corresponding Maxwell-Boltzmann 
distribution, B(E,T) = ^E 2 exp(— E/T)/T 3 , 

(Rppc(E)) = 2^3 J E 2 R ppc (E) exp (-E/T) dE , ( 21 ) 

where the temperature is in units of energy. Combining 
Eqs. (18 1 , (19), and (20 1 with Eq. (21) results in a full three- 


body calculation of the energy-averaged reaction rate, as 
the wave functions in Eq. (20 1 are proper three-body wave 
functions. 

In the special case, where the resonances are very narrow 
and well separated the expression can be simplified greatly. 
If is a bound state, or a very narrow resonance state 
described accurately as a bound state, we can assign a pho¬ 
ton emission width, T 7 , to the transition from this state, 
that is given by [33] 
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where J is the angular momentum of the three-body reso¬ 
nance with energy E R . 

If furthermore the three-body resonance is approximated 
by a Breit-Wigner shape, the photodissociation cross sec¬ 
tion in Eq. (19) can be written as 


<t 7 (.E 7 ) =7t(J + 1/2) 


fhc\ 2 r eff (E)T(E) 

\eJ (. E-E R y + \T{Er 


(23) 




Type is the 


where (T eff (E))~ 1 = (T ppc (E))~ 3 + (T J (E)) 
strong decay width, and the total (in principle energy de¬ 
pendent) width T(E) = T ppc + T 7 . 

With the dominating contributions arising from well sep¬ 
arated, narrow resonances of Breit-Wigner shapes as in 
Eq. (22), the integral in Eq. (21) can be solved analytically, 
and we arrive at a very simple expression for the average 
rate 


(R ppc (E)) = + l lf 2 Veff ^ R) exp (~E r /T). (24) 

{hcphcp, P r /2 t 6 

The only assumption in this expression is that the pho¬ 
todissociation cross section is accurately expressed by the 
Breit-Wigner form in Eq. (22) with three-body resonance 


energy, E R T and width r<I. Other contributing nar¬ 
row well-separated resonances can simply be added. Thus 
Eq. (241 is valid for one or more contributing narrow reso¬ 


nance irrespective of capture mechanism. 

The three-body (two proton plus an even-even core) cal¬ 
culation leads to a 0 + ground state independent of the num¬ 
ber and character of the contributing core-proton single¬ 
particle states. The lowest excited three-body bound or 
resonance state would almost definitely be a 2 + state. This 
is confirmed by the mirror nuclei [55], where the two low¬ 
est excited states are 2 + states. Furthermore, the lowest 
possibly negative parity state is 2.5 MeV above the ground 
state, and there is no 1~ at all in the known spectrum. The 
dominating transition can therefore be assumed to be an £2 
transition. To examine the remote possibility of a 1“ state 
two opposite parity single-particle state must be allowed. 
The relative energies of the 2 + and 1“ states depend on the 
attractions of the two-body potentials for the correspond¬ 
ing partial waves. The energies of these excited states are, 
through Eq. ( [24] ) , all-decisive for the capture rates. Beside 
the three-body resonance energy also the effective width, 
T e ff 7 is important for the capture rate. The realistic as¬ 
sumption that T e f f « T 7 allows an estimate of the relative 
sizes of the 2 + to 1“ effective widths, that is from Ref. [551 


r 7 (g(l+ 1 )) ( 


EryR 


r 7 (£T) \hc{2£+ 3) 


(25) 


where R is the radius of the nucleus. For f=lwe get with 
E 1 = 1 MeV and R = 10 fm that r 7 (£l) « 10 4 r 7 (£2). 
This estimate is from a single-particle model, but within 
the three-body model effects from the other proton can 
at most contribute by a factor of two. The implication is 
that £1 would dominate unless forbidden by rather strict 
conservation laws. 


III. THREE-BODY GROUND STATE 
PROPERTIES 

The nuclear properties of the waiting points are related 
to the borromean structure reached after two-proton cap¬ 
ture. We therefore first investigate the borromean ground 
state depending on the two-body potentials varying along 
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Figure 1. (Color online) The potentials based on the spectra of 
the lowest three A’s for p 3 / 2 (dashed red line) and f 5 / 2 (dotted 
blue line) in isolation, along with the combination of both p 3 / 2 
and f 3 / 2 (solid black line). All two-body potentials have been 
adjusted to produce the energy E 2b = 0.1 MeV. It is seen how 
the spectrum of the combined case always follows the lowest 
available potential. 


the proton dripline. In the following subsections we dis¬ 
cuss three-body energies and structure of the corresponding 
wave functions. 


A. Energies 

The most tightly bound borronrean nucleus is most likely 
a closed (sub)shell for the core while the additional two 
protons occupy empty valence orbits. An extra stability is 
present for N = Z nuclei which conveniently also includes 
waiting point nuclei for N = Z = 32, 34, 36. The ground 
states for both core and borromean nuclei have zero angu¬ 
lar momentum and positive parity as always for even-even 
nuclei. The first excited states of these nuclei are 0.9 MeV 
[35], 0.9 MeV [35] and 0.7 MeV [37], respectively. These 
rather high values demonstrate the stability required for a 
borromean three-body structure build on these inert cores. 

Given the significant Coulomb barrier the three-body 
system does in principle not need to be strictly bound. The 
only requirement is that the proton-decay branch is small 
compared to the /3-decay branch. Small positive three-body 
energies are then possible, which in principle extends our 
energy region of interest. However, the j3 decay half-life 
of 70 Kr has been measured to be 57(21) ms [35], so the 
possible extension is exceedingly small. 

The two-body bound or resonance energies for the dif¬ 
ferent partial waves are now the only pieces of information 
missing before the three-body properties can be calculated. 
We choose to use p 3 / 2 and f 5 / 2 and allow them both simul¬ 
taneously with the same energy as well as one at a time. 
All other partial waves are either not present or shifted to 


Figure 2. (Color online) The three-body energies as a function 
of two-body energies for p 3 / 2 , p 3 / 2 , and / 5 / 2 hi isolation, and for 
P 1/2 and f 3 / 2 as well as p 3 / 2 and f 3 / 2 in combination, the latter 
both with and without p-p interaction, and with and without 
a thr ee-b ody potential. The linear fits are in accordance with 
Eq. (291. All two-body potentials were adjusted to produce 
the same energy, when more that one wave is allowed. Only 
a 0 + state is considered here. The horizontal and vertical line 
indicates our borromean region of interest. 


high energies. 

In Fig. [l] we show the lowest effective adiabatic poten¬ 
tials defined in Eq. (16) for the three cases with the same 
two-body energy, E 2 t, = 0.1 MeV, in all cases. We first 
notice an attractive potential for each contributing partial 
wave. When both p 3 / 2 and f 5 / 2 are allowed we find two at¬ 
tractive low-lying adiabatic potentials. Both contributions 
are present in the combined case where the deepest poten¬ 
tial always follow the lowest potential from either p 3 / 2 or 
fb /2 f° r different p- values. Constructive (or destructive) in¬ 
terference can only marginally occur through non-adiabatic 
terms since p and /-waves cannot couple to form a 0 + state. 
However, this results in avoided crossings for the combined 
case which therefore must produce lower three-body ener¬ 
gies than the individual cases. In general, the combined 
case can then never be less attractive than any of its com¬ 
ponents. 


The coupled set of adiabatic radial potentials is solved 
and energies and wave functions calculated. The two-body 
input parameters are not known although limits from the 
required borromean character can be found. We must 
therefore investigate the three-body properties as function 
of the unknown two-body energies. The three-body ener¬ 
gies, E 3b , are presented in Fig. [2] as functions of two-body 
energies, E 2b , for various selection of contributing partial 
waves. The most spectacular observation is that all curves 
are linear. In the figure we only exhibit results for the most 
interesting energy interval but the same observed linear de¬ 
pendence is accurately followed within the investigated in¬ 
terval, —2.0 MeV < E 2b < 2.0 MeV. 
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When more than one partial wave is allowed, the two- 
body potentials are adjusted to produce the same energy. 
Fixing the energy of one two-body partial wave, while in¬ 
creasing another one from the same value, the E 3b must in¬ 
crease and approach the higher-lying curve corresponding 
to the one contributing wave. If more than two degener¬ 
ate core-proton single-particle levels contribute we would 
find even lower-lying .E^-curves. However, this is highly 
unlikely since this requires nuclear potentials of unprece¬ 
dented high symmetry. Omission of the proton-proton in¬ 
teraction in the case of two contributing partial waves, P 3 / 2 . 
and /5 / 2 , increase the curve even above ( 0 , 0 ), which apart 
from center-of-mass effects should correspond to three non¬ 
interacting systems. Thus, we claim to have established 
limits for the -f^-variation between the lowest curve in 
Fig. [ 2 ] and a parallel curve roughly passing through (0,0). 

So far, the results shown in Fig. [2] have been obtained 
using just the two-body interactions contained in Eqs. ([5]), 
([ 6 ]), and Q. However, it is a well known fact that us¬ 
ing only two-body potentials will usually lead to an un¬ 
derbound three-body structure, as shown for instance for 
17 Ne and 12 C in [23] and [39], respectively. This problem is 
typically overtaken by addition of an effective three-body 
force, V3b(p), to the adiabatic potential given in Eq. (16). 
Nevertheless, inclusion of such three-body force does not 
change the linear dependence shown in the figure. This is 
illustrated by the dashed (yellow) line in Fig. [ 2 } which corre¬ 
sponds to the dot-dashed (black) case, but where a modest 
gaussian three-body force V3 b (p) = S 3 & exp(— (p/po) 2 ) has 
been included. In particular, the values S 3 ;, = — 5.8 MeV 
and po = 6 fm have been used. As seen in the figure, in¬ 
clusion of the three-body force only shifts the three-body 
energy, but keeps the energy relation intact. 


B. Structure 

The simple linear dependency in Fig. [2] is due to the 
special structure of the wave function and the Hamilto¬ 
nian in Eq. <§• The core is much heavier than the proton, 
and the reduced masses are to a very good approximation 
both equal to the proton mass. Then x cp ss y cp and the 
Hamiltonians in Eqs. 0 and 0 are also approximately 
equal, H x ~ Hy. If \^3 b ) denotes the three-body wave 
function, the three-body energy is determined as the ex¬ 
pectation value of the Hamiltonian in Eq. ([ 5 ]), that is 

E 3 b = ('k 3 b | H x |H/ 3 b) + (T 3 (,| H y |H/35) + (^> 3 b\ V pp IH/35) 

— 2 E x + Epp. (26) 

where the approximate equality arises due to the assump¬ 
tion x cp ~ y cp . This approximation is consistent with 
|’I , 3t,) as a product of corresponding two-body wave func¬ 
tions. For instance, when two single-particle states are 
allowed, for example p 3 / 2 and f 5 / 2 , we have Id's;,) = 
A \P 3 / 2 ) x \Ps/ 2 ) y + B | f 5/2 ) x 1 / 5 / 2 ),,, where the angular mo¬ 
mentum coupling to 0+ of both terms implicitly is assumed. 
Product terms of p 3 / 2 and /s/ 2 are not allowed because the 
total angular momentum, 0 + , cannot be reached. 


The above product structure of 1 ^ 35 ) is rather general, 
since the x cp distribution necessarily is described by the al¬ 
lowed single-particle wave functions whereas the same over¬ 
whelmingly dominating single-particle y cp -part for com¬ 
pleteness should be extended to include other angular mo¬ 
mentum components. We then get 

E x = <tf 3 &|fl*|* 36 > = A 2 E 2b (p) + B 2 E 2b (f). (27) 

where the two-body energies are defined by 

E 2 b{p) = (P3/2\H x \p 3 / 2 ) x , 

E 2 b(f) = (f 5/2 \H x \f 5/2 ) x , (28) 

and the cross terms vanish both due to orthogonality of the 
eigenfunctions of H x , and angular momentum conservation 
of the Hamiltonian. We then arrive at the estimate of the 
three-body energy 

E 3b * 2 ( A 2 E 2b (p ) + B 2 E 2b (f)) + <*36|V w |*36> , (29) 


where A 2 + B 2 = 1. The occupation probabilities, A 2 and 
B 2 , are given by the relative weights of the partial waves. 
If only p 3 / 2 or f 5 / 2 are allowed we have A = 1 and B = 0 or 
A = 0 and B = 1, and consequently the linear dependence 
seen in Fig. [2j If E 2b (p) = E 2b (f) the linear dependence 
still results, and in both extremes the slopes of the curves 
are two. Variation between limits must produce a curve 
connecting the corresponding two lines. For example, the 
lower limit, when allowing both p 3 / 2 and /s/ 2 , is given by 
the black, dash-dotted line, and the upper limit is given 
by the curve corresponding to the wave with lowest energy. 
We emphasize the remarkable equality of the slopes of the 
lines in Fig. [2] 

The structure of the solutions is quantified by decom¬ 
position into contributing partial waves. We collect the 
results in Table [T] for the same seven cases as seen in Fig. [2] 
with both sets of Jacobi coordinates. We have chosen two 
rather different values, —2.0 MeV and 2.0 MeV, for the 
two-body energies used to adjust the strength parameters. 
The distributions are very similar for the two energies, and 
the transition from one distribution to the other is slow 
and monotonous. Adding a three-body potential has no 
significant effect on the structure of the wave function. 

The decomposition is rather trivial, when only one single¬ 
particle orbit is allowed. In the second set of Jacobi co¬ 
ordinates, the proton-core set, only this state is allowed. 
However, two couplings to total orbital angular momen¬ 
tum, It = 0 and 1, share the weights, which is equally 
distributed in the first Jacobi set since It is conserved in 
this rotation. 

When two single-particle orbits are allowed the decom¬ 
position includes both structures with different relative 
weights. The It = 0 components are always far larger 
than those of l t = 1. The distribution is strongly depend¬ 
ing on the proton-proton interaction as seen by the com¬ 
plete separation for V pp = 0, where the three-body ground 
state is degenerate corresponding to either |p 3 / 2 ) IP 3 / 2 ) 


° r 1/5/2)^ 1/5/2)^- 

As seen from Eq. (29) the displacement of the lines in 
Fig- 0 not only depends on the proton-proton interaction, 
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Table I. The weights of each partial wave for the same cases as 
shown in Fig. [ 2 ] The state in question is a 0 + state. Here l x 
denotes the relative angular momentum between the two par¬ 
ticles specified by the second column, l y denotes the angular 
momentum of the third particle relative to the center of mass 
of the first two particles, and It is the total angular momentum 
they combine to. The sixth column gives the order, Kmax, of 
the Jacobi polynomium used for the corresponding partial wave. 
The last two columns shows the weight (in percent) of the states 
for different two-body c-p energy. Components where all states 
have a weight less than 10% are generally not included. 








E 2 b 

(MeV) 

Waves 

Jacobi 

lx 

ly 

It K-max 

-2.0 

2.0 

P 3/2 

p-p 

0 

0 

0 

98 

78 

76 



1 

1 

1 

80 

20 

22 


p-c 

1 

1 

0 

80 

79 

77 



1 

1 

1 

80 

20 

23 

P 1/2 

p-p 

0 

0 

0 

98 

48 

45 



1 

1 

1 

80 

49 

52 


p-c 

1 

1 

0 

80 

49 

46 



1 

1 

1 

80 

51 

54 

fs/2 

p-p 

0 

0 

0 

98 

48 

47 



1 

1 

1 

80 

37 

38 



2 

2 

0 

52 

11 

11 


p-c 

3 

3 

0 

74 

60 

59 



3 

3 

1 

74 

40 

41 

P3/2 + /5/2 

p-p 

0 

0 

0 

98 

78 

75 



1 

1 

1 

80 

19 

22 


p-c 

1 

1 

0 

80 

56 

52 



1 

1 

1 

80 

12 

13 



3 

3 

0 

74 

23 

24 



3 

3 

1 

74 

9 

11 

P3/2 + /5/2 

p-p 

0 

0 

0 

98 

80 

77 

(V 36 ^ 0 ) 


1 

1 

1 

80 

18 

20 


p-c 

1 

1 

0 

80 

62 

57 



1 

1 

1 

80 

13 

13 



3 

3 

0 

74 

19 

21 



3 

3 

1 

74 

6 

9 

P3/2 + /5/2 

p-p 

0 

0 

0 

98 

65 

64 

0 


1 

1 

1 

80 

33 

32 


p-c 

1 

1 

0 

80 

66 

65 



1 

1 

1 

80 

34 

33 



3 

3 

0 

74 

0 

0 



3 

3 

1 

74 

0 

0 

Pl/2 + fb/2 

p-p 

0 

0 

0 

98 

66 

61 



1 

1 

1 

80 

31 

35 


p-c 

1 

1 

0 

80 

38 

34 



1 

1 

1 

80 

20 

21 



3 

3 

0 

74 

30 

31 



3 

3 

1 

74 

11 

14 

but also on 

the wave 

function. 

The 

spin- 

■orbit difference 


between pi/ 2 and ps/ 2 causes a slightly different wave func¬ 
tion, which is why the two curves differ. This is also demon¬ 
strated by the fact that displacement changes from —1.54 
to 0.18 MeV, when a very shallow nuclear potential and 
neither spin-spin, spin-orbit, tensor, or Coulomb potentials 
are used in the proton-proton interaction. In fact, by com¬ 
pletely eliminating the proton-proton interaction, the result 
is almost two independent two-body systems. The slight 



CM 



Figure 3. (Color online) The probability distribution for the 0 + 
ground state with equal two-body binding of 0.641 MeV for both 
/ 5/2 and P 3/2 partial waves. Projected contour curves are shown 
at the bottom of each figure. The distance variables correspond 
to the two different Jacobi sets, where the a) and b) panels are 
for the first and second Jacobi sets, respectively. 

positive displacement is in that case due to the fact that 
all three coordinates it, r 2: and r 3 are still coupled in the 
y coordinate as long as the core is not infinitely heavy. 

The actual size of the displacement, caused by the 
proton-proton interaction, can in principle be estimated 
from the matrix element ( v F|V),p|'F). We show the probabil¬ 
ity distribution of the three-body wave function in Fig. [3] 
where the distance of about 5 fm between proton and core 
(Fig. [3 Jd) is a prominent feature. The proton-proton dis¬ 
tance distribution (Fig. [3^,) is much more complicated with 
three peaks at distances of about 2.5 fm, 8 fm, and 11 fm, 
respectively. The Coulomb repulsion would correspond¬ 
ingly be about 0.6 MeV, 0.18 MeV, and 0.13 MeV. 

The strong nucleon-nucleon interaction has strength of 
around 40 MeV and range of 2 fm. To arrive at a total dis¬ 
placement of around —1.7 MeV there must be only a few 
percent of the proton-proton distance-probability within 
the range of 2 fm. To make a reliable estimate of the dis¬ 
placement is therefore very delicate as it depends strongly 
on the solution to the three-body problem. A much better 
computation would be to evaluate E pp directly. However, 
this would still be only an estimate since ’F changes with 
the interaction, and at best we can only reproduce the al¬ 
ready known actual curves in Fig. [2] 
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IV. RADIATIVE CAPTURE 


The critical waiting points in the rp process are defined 
by the long time it takes to capture an additional two pro¬ 
tons in the nucleus. The borromean nature requires a three- 
body reaction producing the strong interaction bound two- 
proton plus core system. In this section we shall focus on 
corresponding reaction rates and the structure of crucial 
intermediate states. The total process is c+p+p —>■ A + 7 , 
which as well can be understood through the reverse pro¬ 
cess, A + 7 —>• c + p +p. It is often very accurate to divide 
part of this process into two steps, A* o A + 7 , where A* 
is one (and sometimes a few) intermediate excited state. 
We shall first investigate the properties of such intermedi¬ 
ate excited states and subsequently calculate the reaction 
rates. Unless otherwise stated a two-body energy of 0.641 
MeV is used, as this is the measured proton separation en¬ 
ergy of 69 Br )5T] . The results are not restricted to such a 
specific energy, but apply rather generally to the region of 
the nuclear chart around the critical waiting points. 


A. Excited continuum states 

The reactions proceed from continuum three-body states, 
that is from two free protons and an (almost) ordinary 
nucleus. This problem can be handled by two conceptu¬ 
ally different methods [39] [42], where the first is to specify 
the boundary conditions directly and solve the Schrodinger 
equation. The second method is to discretize the contin¬ 
uum in a large hyperradial box which in the present case is 
limited by hyperradii less than the box radius, p m ax- 

We shall here use the discretization method with the 
great advantage of using the already defined adiabatic po¬ 
tentials. As discussed in Sec. |II C| the most likely dominat¬ 
ing intermediate angular momenta is 2 + . Also, as discussed 
in Sec. |IIB[ states with the necessary opposite parity to 
form a 1“ state is also very unlikely. However, to later 
study the scale of the El transition, p 3 / 2 and d 5 / 2 waves 
are allowed to form a 1~ state. We show the lowest corre¬ 
sponding adiabatic potentials in Fig.[4]for three such cases, 
where two partial waves for each total angular momentum 
are allowed. The parameters are chosen to be the same as 
already used in the study of ground state properties. 

The potentials in (a) and (b) of Fig. [4] are qualitatively 
similar to those of Fig. [l] with attractive pockets around 
p = 8 fm, infinite repulsion at smaller distances, and bar¬ 
riers at larger p separating regions of interacting and fully 
separated three particles. The attractive region and the 
substantial barrier suggest that there are narrow, low-lying 
resonances or perhaps even bound states. The energies of 
the radial solutions in the box are at first glance also simi¬ 
lar to the ground state 0 + solutions, that is one prominent 
bound, separate, low-lying state and a number of higher- 
lying solutions. However, the bound solutions have posi¬ 
tive energy and would therefore correspond to resonances 
or continuum states. The 1~ potentials shown as (c) in 
Fig- @ are much more repulsive at short distances since 



P [fm] 

Figure 4. (Color online) (a) The potentials based on the spectra 
of the lowest five A’s for the 2 + state allowing both p 3 / 2 and / 5 / 2 
waves. The dashed, horizontal line is at the lowest three-body 
resonance energy. The dotted curve is the sum of Coulomb 
potential and two-body energy, (b) The same for the 1“ state 
with p 3 /2 and d $/2 waves, (c) The same as for (b) with the d 5 / 2 
two-body energy equal to 1.5 MeV. 


the necessary d 5 / 2 state is chosen to be at 1.50 MeV. The 
resonance is at a higher energy, where also more contin¬ 
uum background states have a non-vanishing contribution 
at short distances. 

The three-body energy for the 1~ case with E 2b (p) ^ 
E 2b {d) can be explained using an argument similar to that 
in Sec. ra The three-body wave function can be writ¬ 
ten as 1 ^ 3 b) = A\p 3/2 ) x \d 5/2 ) y + B\d 5/2 ) x \p 3/2 ) y , where 
the angular momentum coupling is to 1~. As m c m p 
the likely weight of the two configurations must be the 
same, and since these are the only possibilities we have 

























Figure 5. (Color online) The three-body energy as a function 
of two-body energy for the 0 + and 2 + state with p 3 / 2 and / 5 / 2 
waves, as well as for the 0 + and 1“ state with p 3 / 2 and d 5 / 2 
waves. The horizontal line is included to guide the eye. 


A 2 = B 2 = 1/2. Then E x = 1/2 (E 2b (p 3/2 ) + E 2b (d 5/2 )), 
and E 3b = E 2b {p 3 / 2 ) + E 2b {d b / 2 ) + E pp . It was found that 
if E 2b (p 3 / 2 ) = E 2b (d 5 / 2 ) = 0.641 MeV then E 3b = 0.61 
MeV, which means E pp ~ — E 2b (p 3 / 2 ) for that specific en¬ 
ergy. However, the three-body wave function is almost un¬ 
changed when changing E 2b (d 5 / 2 ) because both the p- and 
d-wave is needed to form the 1“ state, and E pp is therefore 
also almost unchanged. The result is E 3b ~ E 2b (d 3 / 2 ) for 
E 2 b{P 3 / 2 ) = 0.641 MeV. 

The three-body large-distance configuration is also re¬ 
vealed by close inspection of Fig. [4] A region around 
p = 40 — 50 fm is seen where one potential tend to be 
more flat than dictated by a general Coulomb decreasing 
potential. The size and slow decrease is precisely consis¬ 
tent with the Coulomb potential between one proton and a 
proton-core two-body system in a spatially small resonance 
at 0.641 MeV. This is reflected in the good agreement be¬ 
tween the lowest potential and the corresponding Coulomb 
potential plus two-body energy seen in Fig. [4j In Fig. |4])c) 
the same Vc + E 2b barrier is indicated by the red, dot¬ 
ted curve. This implies that the close-lying proton is in a 
p-orbital. 

The three-body (bound) resonance energies are not nec¬ 
essarily the lowest discretized states. It depends strongly 
on the size of the box since an infinite box must produce a 
continuum of states from the threshold and upwards. We 
detect the resonance state by requiring that the density 
distribution is localized at small distances. The resonance 
energies are shown in Fig. [5] as function of two-body ener¬ 
gies for the choices used for the ground state calculations. 
It is perhaps less surprising to find the same simple and 
accurate linear dependence as we observed for the ground 
states. The criteria for linearity are the same. Also the ef¬ 
fect of a three-body potential would be similar to the effect 


Table II. The weights of each partial wave for the 2 + case, where 
only p 3 /2 and /s/ 2 waves are allowed, and for the 1“ case, where 
only p 3 / 2 and d 5 / 2 waves are allowed. The notation is the same 
as in Table [T] The last three columns shows the weight (in per¬ 
cent) of the lowest resonances. The resonances correspond to 
the first peaks in Fig. [6] Components where all states have a 
weight less than 10% are not included. 


Waves 

Jacobi 

lx 

ly 

It 

Kmax 


Weight 


2+ 

p-p 

0 

2 

2 

199 

27 

8 

15 

7*3/2 + fb/2 


1 

1 

1 

199 

40 

11 

11 



1 

1 

2 

199 

3 

19 

2 



1 

3 

3 

202 

1 

17 

12 



2 

0 

2 

60 

23 

8 

13 



3 

1 

3 

22 

1 

17 

12 


p-c 

1 

1 

1 

160 

41 

12 

10 



1 

1 

2 

160 

38 

4 

1 



1 

3 

2 

162 

3 

15 

15 



1 

3 

3 

162 

1 

20 

23 



3 

1 

2 

162 

3 

15 

15 



3 

1 

3 

162 

1 

20 

23 

2+ 

p-p 

0 

2 

2 

199 

23 

25 

- 

P 3/2 


1 

1 

1 

199 

53 

41 

- 



2 

0 

2 

60 

21 

24 

- 


p-c 

1 

1 

1 

200 

55 

42 

- 



1 

1 

2 

200 

45 

52 

- 

1_ Q 

p-p 

0 

1 

1 

199 

66 

- 

- 

P 3/2 + d 5 / 2 


1 

0 

1 

199 

13 

- 

- 


p-c 

1 

2 

1 

161 

38 

- 

- 

if] 


2 

1 

1 

161 

38 

- 

- 

p-p 

0 

1 

1 

199 

66 

- 

- 

P 3/2 + d 5 /2 


1 

0 

1 

199 

14 

- 

- 


p-c 

1 

2 

1 

161 

38 

- 

- 



2 

1 

1 

161 

38 

- 

- 


a F' 26 (P 3 / 2 ) — ® 2 i,(d 5 / 2 ) — 0.641 MeV. 
b E 2 b (P 372 ) = 0.641 MeV, E 2 b(d b / 2 ) = 1-50 MeV. 


seen in Fig. [2] only it would be more difficult to justify the 
exact size of the potential as the excitation spectrum is less 
well known. 

We can then conclude that the excitation energy is a 
constant independent of the chosen variation in Figs. [5] and 
[2j Specifically we get excitation energies 1.38 MeV and 
1.06 MeV for the 2 + and the 1“ states, respectively. How¬ 
ever, the derived limits of variation are completely differ¬ 
ent. The 2 + excitation energy is strongly limited, since the 
ground state only can move between finite limits as shown 
in Fig. [ 2 ] and the lowest 2 + resonance is entirely determined 
by the p 3 /2 component even when also f b / 2 is allowed with 
the same energy as seen in Table [TT] 

The reason for this behaviour is that any partial wave 
(except angular momentum 1/2) of proton and 0 + -core 
states can be occupied by two protons coupled to both 0 + 
and 2 + . The excitation energy can then vary at most by 
about 1.5 MeV. In contrast, the 1“ state requires both pos¬ 
itive and negative parity proton-core single-particle states. 
Thus, a well defined finite 0 + ground state energy always 
appears whereas the 1“ excitation energy can vary from 
zero to infinity by increasing one of the necessary proton- 
core states towards infinity. 
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The discretization allows normalization of all the states 
in the box. Specifically, as mentioned above the density 
distributions of the resonance states are localized at small 
distances. The probability distributions are almost indis¬ 
tinguishable from that of Fig. [3] each with a peak at a 
proton-core distance of about 6 fm. Only one 1“ reso¬ 
nance is found whereas three rather pronounced low-lying 
2 + resonances appear when both p 3 / 2 and f 5 / 2 are allowed 
with the same energies. By definition all these states have 
the overwhelmingly large probability located at small dis¬ 
tances. The spatial overlaps with corresponding ground 
states are therefore very large in all theses cases. Thus, 
only appropriate angular momentum dependent operators 
and resonance energies are required to initiate highly likely 
transitions, as seen in Eq. (191. 

It is then important to know the angular momentum 
composition of the excited states. The choices of allowed 
partial waves are strongly limiting for these distributions as 
seen in Table [TTJ The simplest are the lowest 2 + state and 
the 1“ resonance. They consists of only proton-core p 3 / 2 
components, and equal measures of proton-core and p 3 / 2 - 
d 5 /2 components, respectively. The two excited 2 + states 
are mixtures of p 3 / 2 and f 3 / 2 proton-core partial wave com¬ 
ponents. 

All discretized continuum states, beside the resonances 
discussed above, are much more dilute and spread out at 
large distances of the box. The spatial overlaps with the 
ground state are therefore very small and any transitions 
would correspondingly be reduced in size. This does not 
necessarily mean that their contributions can be ignored, 
because the number of these states also increase both with 
box size and with energy. At some point they overlap and 
contribute as a genuine continuum. 


B. Cross section and reaction rates 

The three particles in the continuum are not character¬ 
ized by one complete set of quantum numbers. The plane 
wave states for free particles contain all angular momenta in 
a partial wave expansion. In contrast the final nuclear state 
has given angular momentum and parity, and the tran¬ 
sition itself is conveniently specified by a given one-body 
multipole operator. The transitions between well defined 
states are independent of each other, and the prescription 
is therefore to calculate and add the different contributions. 
The transition probabilities decrease strongly with multi¬ 
polarity, which therefore is decided by nature through the 
structure of the borromean final state. 

The 0 + quantum numbers are achieved by coupling of 
the two proton-core angular momenta, which must be un¬ 
occupied by core nucleons. The available low-lying single¬ 
particle orbits therefore depends entirely on the region of 
interest in the nuclear chart. We focuss on the proton 
dripline region around A = 68, and as discussed f 3 / 2 and 
p 3 / 2 are from the mirror nuclei expected to be the dominat¬ 
ing single-particle orbitals, with £2 being the dominating 
transition. For completeness we shall nevertheless investi¬ 
gate the heavily suppressed El transition. 



Figure 6. (Color online) The cross sections for 2 + —► 0 + with 
both p 3 /2 and / 5 / 2 waves (solid, black curve), and only with p 3 / 2 
waves (dashed, blue curve). Also included are the cross section 
for 1" ->■ 0+ with P 3 /2 and d 5 j 2 - For the red dotted curve 
Ed 5/2 = 0.641 Mev, and Ed 5/2 = 1.50 MeV for the green dash- 
dotted curve. In both cases the El transtition has been scaled 
down by 10 3 to make the figure readable. The photodissociation 
energy is related to the total energy, E, and the binding energy, 
B,by E^ = E+\B\. 


We proceed by calculating the discretized three-body 
continuum states for given total angular momentum and 
parity with selected two-body input energies. The cross 
sections and reaction rates can then be obtained by summa¬ 
tion over these discrete “continuum” states. The method 
does not assume any specific reaction mechanism, that is 
two-step via a photon emitting resonance and/or contin¬ 
uum state, and both sequential and direct reactions are 
included in the numerically obtained results. 


The discretization implies that the low-energy spectrum 
can be too sparse when the box size is comparable to the 
extension of the potential barriers. The space outside the 
barriers is then too small to provide box bound states. On 
the other hand, an attractive short-distance pocket would 
produce an isolated bound-state like resonance which then 
would mediate all the low-energy transition probability. 
The missing continuum states in its energy neighbourhood 
would have vanishing spatial overlap with the ground state, 
and consequently also vanishing transition probability. In 
the high-energy limit the level density increases and the 
contributions are distributed over many levels. At some 
point the included Hilbert space in the basis becomes insuf¬ 
ficient. Fortunately, we are able to cover an energy region 
sufficient for the astrophysical reactions of interest. 

In Fig. [6] we show the cross section from Eq. (19) as 
function of energy above threshold for two cases of 2^ and 
1“ excitations each. All the peaks are at the resonance 
energies where the spatial overlaps to the ground state allow 
finite cross sections. Two clear differences are seen between 
the two 2+ cases. First of all, the two lowest p 3 / 2 + f 5 / 2 
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peaks are pushed a little bit further apart as for a two- 
level system with an additional interaction. Secondly, the 
peaks are notably smaller when only p 3 / 2 is allowed. The 
reason is the partial wave composition of the peaks and the 
number of contributing potentials. 

All peaks are in the £> 3/2 case composed entirely of 
( l x ,l y ) = (1,1), as seen in Table |H} where only one po¬ 
tential contributes significantly to the lowest peak. In the 
P 3/2 + fd /2 case the same structure is found in the low¬ 
est peak where two potentials now contribute evenly. The 
next two peaks are in contrast composed of an almost even 
mixture of (l x ,l y ) = (1,3) and (3,1), which clearly is not 
allowed with only contribution from the p 3 / 2 state. In be¬ 
tween the resonance peaks are the part of the cross section 
which makes the non-resonant contribution to the reaction 
rate. This contribution is completely negligible compared 
to the resonant contribution for low energies while increas¬ 
ingly appearing at higher energies. As the lowest peak con¬ 
sists almost exclusively of p-waves, adding higher angular 
momentum orbitals would not change the cross section nor 
the resulting reaction rates. 

The cross section for the 1“ excitation is also shown in 
Fig. [ 6 ] When the d 5 / 2 and p 3 / 2 two-body states have the 
same energy we find only one huge contribution at very 
low energy corresponding to one and only one resonance. 
When the d 3 / 2 energy is increased this cross section peak 
(resonance energy) moves up in energy, gets broader, and 
decreases in size. It is found that increasing the d 3 / 2 en¬ 
ergy to 3.00 MeV makes the 8 1 contributions insignificantly 
small compared to the 82 contributions in Fig. [ 6 ] We can 
conclude that 8 1-transitions are large even for relatively 
high-lying two-body states which are necessary for the com¬ 
position of the 1 “ continuum or possible 1 “ resonances. 

The capture process takes place in an environment where 
temperature is an important parameter. The three-body 
energy is therefore not a priory given, but occurs with a 
certain probability distribution and with the capture rate 
specified in Eq. (21 1 . The resulting rates are given by the 
full lines in Fig. Jj The lowest resonance peaks in Fig. [ 6 ] 
suggest narrow states, which can be approximated by a 
Breit-Wigner shape. It is then possible to use the much 
simpler expression in Eq. ( |24| ) to find the contribution from 
each resonance to the overall reaction rate. 

Figure [7] shows that the lowest resonance clearly dom¬ 
inates for both cases even at temperatures well above 3 
GK. Summing the contribution from the isolated lowest 
resonances results in a better agreement with the full cal¬ 
culation to a higher temperature, but large deviations ap¬ 
pear above around 5 GK. It is very remarkable that such 
a simple expression, based on a single resonance, is able to 
estimate the reaction rate so accurately over several GK. 
It is even more impressive that a simple sum over the res¬ 
onances can further increase the temperature range of its 
applicability while maintaining the accuracy. Increasing 
the g? 5/ 2 energy increases the resonance energy, which in 
turn increases the temperature with highest rate, while si¬ 
multaneously reducing this rate. The characteristic peak 
in the reaction curve is also smeared out by the effects of 
higher energy continuum contributions. 



Temperature [GK] 


Figure 7. (Color online) The reaction rates corresponding to 
the cross sections from Fig. [6] The result of the full calculations 
is given by the solid curves. The dashed curves in the same 
colour are the result of applying Eq. (241 to the lowest peak 
in Fig. [6] The dotted curves are the sum of contributions for 
the three (two) lowest resonances for p 3 / 2 + /s /2 (P 3 / 2 ) using 
Eq. ( |24| ) . The 81 rates have been scaled down by 10 4 to make 
the figure readable. The various resonance widths for the lowest 
resonances can be found in Table Hill 


There are a number of possible corrections, which should 
be considered. First of all, the rates in Fig. [7] are based on 
the full three-body 2 + ( 1 ~) spectrum in the energy region 
seen in Fig. [bj If, as in the mirror nuclei (34] . several of the 
lowest states are 2 + states, they are all included. Contri¬ 
butions from other possible states, such as 4 + or 3“ states, 
could conceivably contribute, but they would be suppressed 
by several orders of magnitude because of the higher order 
transition necessary. All transitions are to the same 0 + 
ground state. Excited, bound three-body states are not 
accounted for as they are very unlikely at the edge of the 
dripline, in particular considering the lowest excited state 
in the mirror nuclei is 0.94 MeV above ground (Ml . 


Due to the high temperatures involved core excitations 
could also be a contributing factor. The first excited state 
in 68 Se is E cl = 0.854 MeV above the ground state [351 . 
The probability of occupation at a temperature of 4 GK 
would then be exp (— E c \/(ksT)) =0.09. To get the reac¬ 
tion rate from Eq. ( | 21 [ ) we must multiply by R cpp which 
is proportional to the photodissociation cross section as 
parametrized in Eq. (23) for a resonance. The peak struc¬ 
ture suggested by the parametrization is used to derive 
Eq. (24) which assumes that the barrier penetration fac¬ 
tor contained in T and T e f f has a relatively smooth energy 
dependence. The numerical calculations of <r 7 shown in 
Fig. [ 6 ] reveal the expected peak structure for small E 7 
confirming the sufficient smoothness of T and F e f f. With 
everything being the same for ground and excited states 
we conclude that the contribution from low-lying excited 
states are reduced by their Boltzmann factor which then 
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would be only a minor contribution. The second excited 
core state is 1.20 MeV above ground, and is even less sig¬ 
nificant. To include contributions from core excitations or 
other bound three-body states one would have to include 
the partition functions of both the initial target and final 
three-body nucleus 01m, to account for the thermal equi¬ 
librium between the ground and the excited states. This 
would in this case be minor corrections. 

The discretized three-body states are limited to energies 
less than 10 MeV, so the rates must decrease for sufficiently 
high temperatures. Going to higher temperatures the com¬ 
puted £2 transition rate has a maximum value of about 
6.5 and 3.5 x 10~ 11 N'^ cm 6 s“ 1 mol _2 at just above 25 GK 
for the black and the blue curve respectively, and then de¬ 
creases similarly to the £ 1 transition. The maximum value 
for the £1 transition occurs at much lower temperatures 
as there, as seen from Fig. [6j is no significant contribution 
to the cross section at higher energies. In both cases, the 
restriction imposed by the limit in discretized energies is 
only relevant far outside our region of interest. 

As the single peak approximation fails at higher tempera¬ 
tures it would also fall short at lower temperatures. At low 
temperatures even the lowest resonance state would not 
have a significant probability of being populated, and the 
rate would be dominated by off-resonance contributions. 
Fortunately, this is well below our region of interest. 

If both a 2 + and a 1“ state was present simultaneously 
transitions like 2 + —► 1“ —> 0 + or 1“ —> 2 + —> 0 + could be 
imagined, which would affect the final reaction rate. This 
effect could be computed by establishing the relations be¬ 
tween the relevant two- and three-body energies, and then 
calculating the reaction rate in the same way as before. As 
before, the effect of any higher lying resonances would only 
be significant if the temperature was sufficiently high, as 
demonstrated by the accuracy of the lowest peak in isola¬ 
tion in Fig. [7] For our purposes it is therefore unnecessary 
to consider such corrections. 

If a three-body force is added and the energy levels 
shifted, the effect on the cross section and reaction rate is 
only an energy scaling determined by Eqs. (19) and (24) re¬ 
spectively. The overlap matrix element would not change, 
as the energy shift does not change the structure of the 
wave function as seen from Table U 


V. PRACTICAL IMPLICATIONS 


The final purpose of this paper is to provide a method 
for estimating proton capture reaction rates for the region 
around A ~ 70 and N ~ Z in general, and the critical wait¬ 
ing points in particular. To that end we first discuss the 
reaction mechanism resulting from the sets of input param¬ 
eters. Second, we present estimates of the crucial two-and 
three-body ground and first excited state energies, and fi¬ 
nally, we combine these considerations to present a general 
method for estimating radiative capture rates around the 
critical waiting points. 


A. Reaction mechanism 

The calculations are carried out without need for spec¬ 
ification of the reaction mechanism. In contrast, from the 
calculated results we can deduce how the process proceeds 
between the borromean bound state and the three free 
constituents in the continuum. Due to the principle of 
detailed balance both reaction directions are equally well 
suited for both qualitative and quantitative descriptions. In 
words the process starts with bombarding the borromean 
bound state with photons of energy larger than the three- 
body binding energy, where the populated continuum states 
somehow end up as three independent infinitely separated 
three particles, that is two protons and a core. 

The process A + 7 — > p + c + p could proceed in several 
possible steps. It could include an excited continuum state 
of the borromean nucleus, A*, or a quasi-stable two-body 
proton-core (pc) configuration, although either or both of 
these steps might be skipped. In this way the process is di¬ 
vided in distinct, and possibly, independent steps. Whether 
they are followed or not in a sequential progression is defin¬ 
ing the reaction mechanism, and is used to categorize it. 

The traditionally denoted sequential path [43} [44] is with 
the proton-core state, but with or without A*, 

A + 7 —» A* —►( pc) + p —> p + c + p, 

A + 7 —► (pc) + p^p + c + p. 

Likewise, the route usually called direct mm, again with 
or without A*, is 


A + 7 —1 A* —>p + c + p, 
A + 7 —> p + c + p. 


The formulations and calculations presented in the pre¬ 
vious sections are schematic or quantitatively accurate de¬ 
pending on the point of view. The schematic impression 
arises from the relatively strong assumptions of the few 
contributing proton-core partial waves. This is, however, 
not unrealistic within the field of few-body cluster models, 
which has also been able to provide an appropriate descrip¬ 
tion of two weakly bound nucleons surrounding an ordinary 
nuclear core. The quantitative accuracy emerges as soon as 
the few-body model approximations are shown to be cor¬ 
rect, because then the reliability and completeness of the 
results are unavoidable. 


The path chosen by nature depends on the characteristics 
of the system. However, substantial numerical simplifica¬ 
tions as well as insight can be gained with a dominating 
narrow resonance. Approximating the cross section by a 
peaked function like a Breit-Wigner shape we arrive at 
the extremely simple temperature dependent rate expres¬ 
sion in Eq. (24). This is exactly the same expression as in 
Eq. (15) of Ref. (32], with F 7 replaced by £ e ff- However, 
in Ref. |32] this reaction rate is derived in the “extreme 
sequential limit”, under the much stronger conditions that 
r 7 is much smaller than T ppc , and the direct decay is dis¬ 
allowed. If T 7 -C T ppc then T eff ~ r 7 , and the expres¬ 
sions become identical. On the other hand, if T ppc -C F 7 


12 





f N 

"S' 

Q. 



‘l/l 


12 

9 

6 

3 

0 



12 

10 

8 

6 

4 

2 

0 


Figure 8. (C olor online) The square of the angular wavefunction 
<f> from Eq. (121, multiplied by the phase factor cos 2 (a)sin 2 (a) 
and integrated over and Q y , as a function of p and a for the 
lowest A of the 2 + state in the second set of Jacobi coordinates. 


then r e ff — T ppci and photon emission dominates over the 
strong decay channel. 

Independent of validity of the simplified rate expression 
in Eq. (24), the process can still be either direct or se¬ 
quential, or for that matter any mixture. This is perhaps 
better appreciated by explicitly explaining that T 7 is de¬ 
termined entirely by the excited and ground state short- 
distance properties, while the strong decay is entirely de¬ 
termined by the structure underlying the barriers which 
has to be overcome before the reaction is completed. The 
effective thickness of the barrier is determined by the three- 
body resonance energy level. 

The intermediate structure can be directly investigated 
by the density distribution of the corresponding angular 
wave function from Eq. ©• The angular wave function 
corresponding to the lowest A n is shown in Fig. [8] as a func¬ 
tion of p and a pc . The structure is abundantly clear for 
p values larger than about 20 fin. The probability distri¬ 
bution is only finite at large and small a pc , which through 
Eq. 0 implies that the distance is small between either 
the core and the first proton, or the second proton and the 
centre of mass of the core-proton system. This is precisely 
the proton-core resonance configuration properly antisym¬ 
metrized. As this is the configuration for the lowest A n , the 
energetically most advantageous escape route is apparently 
through this configuration corresponding to sequential de¬ 
cay. However, this conclusion may change with relative 
sizes of the two- and three-body resonance energies, and as 
function of the total three-body energy. The angular wave 
functions corresponding to higher lying A„ eigenvalues have 
different configurations. By including all the relevant, rela¬ 
tively low-lying A spectrum, all the relevant angular config¬ 
urations, and thereby all the relevant reaction mechanisms, 
are included, and not just the extreme sequential or direct 
paths. 

The decay width from resonance to three particles can 
be estimated by use of the WKB tunnelling probability 
through the lowest potential barrier. This is most accu¬ 


Table III. The three-body resonance energy along with both 
photon decay width calculated using Eq. (22 ) and proton decay 
width based on a WKB calcualtion. The effective width F e f f is 
calcualted from these values. The cases of interest correspond 
to the peaks in Fig. [6] All energies and widths are in MeV. 


State 

Er 


r 

ppC 


r 7 


be// 


2 + : p j + /j 

1.14 

4.1 • 

10 _a 

5.3- 

10“ 1U 

4.7- 

10=™ 


2.17 

5.9- 

IO -3 

8.3- 

IO' 10 

8.3- 

lO-w 


4.29 

3.0- 

10 1 

2.7- 

10 _s 

2.7- 

IO -8 

2+ : p-i 

1.38 

2.4 • 

io - 7 

5.4- 

io- 10 

5.4- 

IO” 10 


3.68 


9.6- 

lO -1 

1.5- 

10" 8 

1.5- 

10" 8 

1“ : pa + dj 

0.61 

a 



1.7- 

IO" 6 


- 


1.5C 

1> 

5.2- 

IO" 7 

1.0- 

IO -5 

4.9- 

10" 7 


a ®26(P3/2) — E 2 b(d$/ 2 ) — 0.641 MeV. 
b E 2b (p 3/2 ) = 0.641 MeV, E 2b (d 5/2 ) = 1-50 MeV. 


rately done using numerical integration between the classi¬ 
cal turning points where the effective potential equals the 
three-body resonance energy, V ef f(p 0 ) = V e ff(p t ) = E R . 
Based on the energy difference from the potential minima 
to the resonance energy, the harmonic oscillator frequency 
uj o can be determined. The knocking rate from Ref. [33] is 
about lo/tt which must be multiplied by the second order 
WKB tunnelling probability, 1/(1 + exp(2S)), that is 


Tppc = —- (1 + exp(2S)) 1 

7r 


(30) 


5=1 


rPt 


\J ‘ZlApcp ( VeffiP ) - E R ) d P ■ 


The T 7 width can be calculated using Eq. (22). The 


effective width T e ff can be calculated from T 7 and r ppc . 
The values of these key quantities for the three-body de¬ 
cay, relating to the cases of practical interest studied in the 
previous sections, are presented in Table |III| We find in 
almost all our borromean cases that A* is a well defined 
resonance with a very narrow width. The reason is readily 
found in the potentials, where the attraction is sufficiently 
strong, inside very pronounced confining barriers, to hold 
a narrow resonance. As seen in Eq. (22) T 7 depends on the 
resonance energy, but it is dominatea by the overlap ma¬ 
trix element for the given multipole transition. As such it 
increases slowly with energy irrespective of partial waves, 
but it changes dramatically when the order of the transi¬ 
tion is changed. The increase is much more drastic for T ppc 
as it depends exponentially on barrier thickness. 

When the energy of the 1“ state is increased, by increas¬ 
ing the energy of the d 5 / 2 single-particle proton levels, the 
resonance becomes wider. The attractive pockets of the 
potentials diminish while moving to higher energies, imply¬ 
ing that the resonance features disappear. The continuum 
states would then have substantially less relative probabil¬ 
ity at distances comparable to the size of the ground state. 
This decreases the rate of the 1“ —> 0 + transition. We 
emphasize that the method of the numerical calculations is 
completely unchanged, and the results are in fact obtained 
without any information about whether such a resonance 
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exists or not. 

For the lowest resonance the r ppc is only about an order 
of magnitude larger than T 7 . This implies the lowest res¬ 
onance is at the edge where the effective width becomes a 
mixture of the two. As T ppc decreases exponentially with 
En without lower limit, while T 7 decreases as E %~' +1 , where 
the photon energy necessarily is finite, F e // will be domi¬ 
nated by T ppc if the resonance energy is lowered. 

If the three-body resonances, Eh, are the only interme¬ 
diate configurations acting as doorway states, then non¬ 
vanishing processes occur only when E ss Eh where the 
uncertainty is determined by the width of the states. On 
the other hand, small rates and cross sections arise also for 
energies without the match to resonance doorway states. 
The rate is calculated from the theoretical formulation for 
given energy, but the applications are rather for given tem¬ 
perature, T. The temperature smearing over the energies 
up to around T then all contribute, and the higher energy 
contributions become exponentially suppressed. A much 
more detailed discussion can be found in Ref. [32.. 

The character of the reaction mechanism as sequential, 
direct, a mixture or something else, is fundamentally de¬ 
termined by the dynamic evolution from small to large 
distances. However, a number of rather solid conjectures 
can be made from the calculated static properties. This 
can be elucidated with two schematic potentials based on 
extreme geometric progression corresponding to sequential 
and direct decay, for more details see Ref. [35] . Outside 
the strong attractive region the Coulomb potential for one 
proton moving away from proton-core state of given pos¬ 
itive energy, E pc , is given by V seq = (Z c + l)e 2 /p + E pc . 
In case of several possible, different two-body energies the 
most likely potential to tunnel through would be the small 
and narrow potential. If the two protons are moving away 
from the core along the most favoured symmetric linear 
configuration with the same p the potential is instead, 
Vdir = (2 Z c + l/2)e 2 v / 2/P, where Z c |e| is the core charge. 
These two potentials, V seq and Vdir, cross each other at p c 
for an energy V c , that is 


e 2 ((2i/2 — 1)Z C — 1 + 1 /a/2) 
— e 

^pc 

v _ p (2Z C + 1/2)V2 

c PC (2V2-1)Z C - 1 + 1/72 


/ 2.6Z C 
E pc 

1.6E pc 


(31) 

(32) 


for Z c 1. With E pc = 0.641 MeV and Z c = 34, we get 
p c = 139 fm and V c = 1.00 MeV, which is significantly fur¬ 
ther out than the potential thickness in Fig. [4jr. The most 
energetically favourable decay path is then the sequential 
path, barring drastic, and most likely energy demanding, 
changes in the configuration during the decay process. As 
Em ~ E 2 b for l - the potential thickness is tending towards 
infinity, which would exclude a sequential decay. 

The reaction mechanism for photo dissociation, or equiv¬ 
alently radiative capture, of given energy E can then be 
expected characterized as sequential for V c < E, mixed se¬ 
quential direct for E pc < E < V c , mixed direct and virtual 
sequential E < E pc < V c , and direct E « E pc [44l|46j. We 


emphasize that these characteristics are not rigorous prop¬ 
erties, although limiting cases would be observable m , as 
they would leave distinct signatures in the energy distribu¬ 
tion of the emerging fragments. Based on these limits the 
reaction mechanism can very convincingly be classified as 
sequential for the 2 + case, direct for the low-lying 1” case, 
and either sequential or a mixture for the high-lying 1“ 
case. It should be noted that this is very strongly depen¬ 
dent on the specific potential depths and the temperature 
in question. For temperatures significantly higher or lower 
than the resonance energy the reaction mechanism would 
not be determined by the single lowest resonance, but by 
the continuum background contribution or possibly a com¬ 
plicated combination of several different resonances. 


B. Energy level predictions 

The central parameters in the effective lifetimes of the 
critical waiting points are the proton binding energies of 
the two following isotones. Ideally, one would like to predict 
these two- and three-body energies exactly. Unfortunately, 
insufficient experimental knowledge makes such predictions 
difficult. However, as seen in Sec. |III A[ it is possible to 
establish relations between the possible two-body energies 
and the needed three-body energies. From there limits can 
be inferred regarding the position of low lying excited en¬ 
ergy levels in both the two- and three-body system. 

As the critical waiting points are borromean in nature 
the core-proton system must be unbound, while the two- 
proton system must be bound. In other words, the region 
of interest is limited to E 2 b > 0 and E 3 b < 0. The following 
limits are calculated without including a three-body poten¬ 
tial. Adding an attractive three-body potential would lower 
the three-body energy, and thereby increase the estimated 
energy ranges. 

If the ground state proton separation energy, S p , is 
greater than 0.35 MeV (where the highest single wave curve 
in Fig. [2]crosses zero), then there must be a very close lying 
first excited core-proton state for the three-body system to 
be bound. Depending on the S p value a very narrow energy 
range for this excited level can be predicted. Likewise, S p 
cannot be greater than 0.74 MeV (where the p 3 / 2 + / 5 / 2 
curve crosses zero in Fig. [2j) , as the three-body system then 
could not be bound. Of course, this is based on the as¬ 
sumption that shell model predictions and the mirror nu¬ 
clei correctly identifies the relevant single-particle states in 
the given region. 

The only thing missing, if proper estimates are to be 
made concerning the waiting point nuclei, is the two-body 
energies. However, even this is not as severe a limitation as 
might be imagined. Based on the results in Fig.[2]energy in¬ 
tervals can be established for the proton separation energy. 
Recently, it was possible to measure the ground state pro¬ 
ton separation energy of 69 Br to S p ( 69 Br) = 641(42) keV 
[4l] . For this to comply with our result the ground state 
must be either a p 3 / 2 or a / 5 / 2 state, with the other be¬ 
ing a low lying excited state, otherwise the three-body sys¬ 
tem would be unbound. In Ref. [41] the ground state is 
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surmised to be a f 5 / 2 state with a p 3 / 2 state lying an un¬ 
known distance above. All lines allowing only one partial 
wave have crossed into the unbound three-body region be¬ 
fore E 2b = 0.641 keV. The same is true for the f 5 / 2 +P1/2 
line. The only remaining possibility among the likely par¬ 
tial waves is the p 3 / 2 + f§/ 2 combination. 

The upper and lower limits for this first excited state 
can now be established from the relation between two- and 
three-body energy. The lower limit is of course E 2 b{p 3 / 2 ) = 
0.641 MeV, corresponding to a degenerate ground state. 
The upper limit is determined by keeping the f 5 / 2 energy 
constant and varying the p 3 / 2 energy. Setting E 2 b(f$/ 2 ) = 
E 2 b{p 3 / 2 ) = 0.641 MeV results in E 3 b = —0.21 MeV. By 
slowly changing E 2 b(P3/ 2 ) the upper limit is found to be 
at E 2 b{p 3 / 2 ) = 0.800 MeV, where E 3 b = —0.00 MeV. The 
first excited p 3 / 2 state in 69 Br is then predicted to lie less 
than 0.16 MeV above the ground state. Likewise, using the 
estimates included in the AME 2012 collection, a proton 
separation energy of 5 p ( 73 Rb) = 0.6 MeV is predicted [45] . 
This implies the ground and first excited states consist of a 
P3/2 and fb /2 state, one lying no more than 0.2 MeV from 
the other. For 65 As the AME estimated proton separation 
energy is at 0.09 MeV, which makes it difficult to predict 
anything specifically. However, it does make it very unlikely 
that there is a p\/ 2 or a f 5 / 2 ground state without a low 
lying first excited state. 

Similar predictions can be made concerning the excited 
continuum levels. Based on Fig. [ 5 ] the 2 + state is found to 
very consistently be 1.38 MeV above the ground state. The 
lower limit of the resonance energy, Er, then corresponds 
to the lower limit of the ground state, i.e. = 2 x 

0.641 MeV - 1.54 MeV + 1.38 MeV = 1.12 MeV. Likewise, 
the upper limit corresponds to where the ground state is 
at the edge of being unbound, i.e. E 2 b = 0.74 MeV and 
Er &x = 1.38 MeV. This very narrow interval can be used 
to estimate limits for the reaction rate. 

It is considered much more unlikely that the three-body 
system will form a L state based on both mean field cal¬ 
culations and comparisons with mirror nuclei. However, as 
has been shown, if it is even remotely possible the dipole 
transition will dominate the reaction rate. The placement 
of the 1“ resonance is not as sharply limited as the 2 + reso¬ 
nance, which is dictated by the lowest two-body resonance. 
On the other hand, the 1“ resonance needs a combination 
of opposite parity states, and could move to arbitrarily high 
energies without affecting the 0 + ground state. 


C. Rate predictions 


In the preceding sections it was argued that very narrow 
resonances will be produced by the fairly attractive short- 
distance region in combination with the wide Coulomb bar¬ 
rier at large distances. In addition, very confining lim¬ 
its have been placed on both two- and three-body energy 
levels. Collectively, this allows for limits to be placed on 
the relevant three-body reaction rates for specific energies. 
These limits would again be affected by the addition of a 


three-body potential. However, the energy scaling is pre¬ 
dicted by Eqs. (19) and (24). More generally, this provides 
a method for estimating proton capture rates for the region 
in general given few experimental data. 

As the reaction rate is dominated by narrow resonances 
it can well be approximated by Eq. (24) for temperatures 
in the 0.1 — 4 GK range. The limits of the resonance en¬ 
ergy was established in Sec. [VB] for the specific case of the 
critical waiting points. Assuming the overlap matrix ele¬ 
ment is constant for small changes in two-body energy, T 7 
only depends on photon energy as E^ l+1 . However, as the 
distance between the curves in Fig. [5]is constant, E 1 would 
also be constant independent of two-body energy. The same 


can 


is therefore true of r~, and the T 7 value from Table III 
be used for a given resonance. As T= r ppc + r~ 1 the 
largest possible value of T e // is the constant F 7 . 

It is then clear that using Eq. (24) the largest rate is 
achieved at the lowest resonance energy where F e ff ~ T 7 . 
At some point, for sufficiently low resonance energies, the 
r e // = r 7 assumption no longer holds as r ppc continues to 
decrease. For lower resonance energies F e ff < T 7 and the 
rate decreases. 

To estimate the rate it is then necessary to estimate r ppc 
which means the barrier must be determined very accu¬ 
rately. Unfortunately, the intermediate distances that are 
relevant here are notoriously difficult to treat accurately 
by simple expressions. For large distances approximation 
such as the V,n r and V seq potentials presented earlier could 
be used. These will however overestimate the potential, 
which will underestimate of the rate exponentially. A more 
appealing alternative is to use the potential calculated in 
full in Sec. |IV A| The potential outside the barrier is the 
determining part. To study the rate based on a two-body 
energy different from the value of 0.641 MeV this part of 
the potential needs to be shifted an amount corresponding 
to the difference in two-body energies. 

The final procedure in estimating the three-body reaction 
rate is as follows. The first step is to asses the two-body 
energies somehow. For waiting point nuclei limits can be 
established from Fig. [2] The second step is to find the 
three-body energies based on the two-body energies. The 
ground state three-body energy is determined from Fig. [2j 
while the first excited level is determined from Fig. [5] If 
different two-body energies are used for the relevant partial 
waves the lower limit for the three-body energy is given by 
the curve where the energy is the same, while the upper 
limit is given by the highest single wave curve. The third 
step is to estimate T e f f. The value of T 7 is given in Table 
m for the most likely cases. Otherwise the single-particle 
decay rates can be used to estimate T 7 [33], where a factor 
of two should be added to account for the two particles. In 
many cases T 7 <C r ppc , and the exact value of r ppc is then 
not relevant. When necessary, for relatively low resonance 
energies, r ppc can be estimated with a WKB calculation 
through the known potential barrier shifted appropriately 
according to the two-body energy. The fourth and final step 
is to calculate the limits of the reaction rate with Eq. ( [24] ) 
using the calculated effective width and the limits on the 
resonance energy. This method applies generally for proton 
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Figure 9. (Color online) Estimates of the reaction rate for 2 + —► 
0 + with p 3 /2 and / 5 / 2 based only on two-body energies. The 
blue curve at E 2 b = 0.74 MeV corresponds to the upper limit 
established in Sec. |V B| The low temperature region has been 
scaled up to demonstrate the similarities. 


capture in the region around the critical waiting points. If 
considering 64 Ge or '°Kr instead of 68 Se minor changes in 
the long-range Coulomb potential would have to be consid¬ 
ered, but the most important change would be the two- or 
three-body energy spectrum. The same method could be 
applied based on estimates of the particular energy levels. 

The result is seen in Fig. [9J where the 2 + —► 0 + rate is 
estimated based on three different two-body energies. The 
energy studied in Sec. |IV A| is used along with the upper 
energy level of E 2 b = 0.74 MeV established in Sec. VB and 
an energy an equal amount below. The low temperature 
region is scaled up to show how closely the rates agree for 
temperatures below ~ 1 GK. The rate decreases with in¬ 
creasing resonance energy because of the exponential factor 
in Eq. (|24|) , and it also decreases with decreasing resonance 


energy because we are at the edge where T ppc starts affect¬ 
ing T ef f. Changing the various two-body energies individ¬ 
ually, between 0.54 — 0.74 MeV, would change the three- 
body resonance energy, but the final reaction rate would be 
within the limits here presented. Going outside this two- 
body energy region would decrease the rate exponentially, 
based on Eq. (24), for temperatures around 1-4 GK. 

These limits are based on the assumption that only the 
lowest resonance contributes significantly. The tempera¬ 
ture must be in the vicinity of the three-body energy for 
this to be true. For much lower temperatures the reso¬ 
nance energy is not accessible, and the main contribution is 
through off-resonant reactions. This would lower the rate 
exponentially. For higher temperatures other resonances 
would contribute, and the lowest resonance might only con¬ 
tribute through the tail of the cross section. 

In principle, these rate estimates should include the var¬ 
ious contributions discussed at the end of Sec. IV B[ even 


for temperatures in the 1-4 GK region. However, these are 
all minor corrections which does not change the result sig¬ 
nificantly. Major changes would only occur with changes 
in energy levels. 


VI. CONCLUSION 

The waiting points in the rp process are prime candi¬ 
dates for three-body calculations by their very definition. 
Given the heavy core in the three-body system the proton- 
core interaction will determine the characteristic of the 
system, with only minor corrections stemming from the 
proton-proton interaction. Because of this, very simple re¬ 
lations between the two- and three-body energies can be 
derived and understood based on the three-body Hamilto¬ 
nian. These relations apply equally well for both ground 
and excited energy levels. These very simple, and rather 
general relations between two- and three-body structures 
and energies allow for the estimation of either two- or three- 
body properties given very sparse experimental data. 

It is also seen that the three-body structure consists of 
the two protons located at the surface of the heavy core. 
At larger distances the structure is that of single proton 
moving away from a two-body core-proton resonance. Con¬ 
ceptually, this corresponds to what is traditionally known 
as sequential decay. For low available initial three-body 
energy (and low temperature) this configuration would no 
longer be the energetically most favourable. 

When calculating the cross section for the process A + 
7 —> p + p + c the effective potential confining the parti¬ 
cles is decisive. The relatively heavy core gives rise to a 
large Coulomb barrier and an attractive short-distance po¬ 
tential. This makes the lowest continuum states narrow, 
well-defined resonances. Based on the spectrum of the mir¬ 
ror nuclei, 2 + states were assumed to dominate the low- 
lying spectrum. The full three-body rate calculation was 
based on the entire 2 + spectrum. Assuming the narrow, 
well-defined resonances have a Breit-Wigner shape, the full 
rate calculation could be reproduced from around 0.5 GK 
and up to 4 GK based solely on the lowest resonance. Sum¬ 
ming over the contributions from the first few resonances 
for the 2 + case increased this temperature range to about 
5 — 6 GK. The off-resonance, background contribution is 
mainly relevant for lower temperatures, where the chance 
of accessing the resonance level is much smaller. For the 
much more unlikely 1“ case only one resonance state is 
present, which then dominates a wide temperature range 
around the resonance energy. 

In the unlikely scenario where an 1“ state is available, de¬ 
pending on its position, the £1 transition could contribute 
significantly. However, based on both mean field calcula¬ 
tions and experimental measurements in the region around 
the critical waiting points single-particle orbitals of like par¬ 
ity dominate the low-lying energy levels. Also no 1“ state 
is seen in the mirror nuclei. This very effectively excludes 
three-body states of negative parity. The most likely tran¬ 
sition is then the always allowed (as long as more than 
s-waves are available) £2 transition. 


16 























As the expression based on the assumed Breit-Wigner 
shape depends mainly on the resonance energies, the re¬ 
lations between two- and three-body energies can be used 
to provide limits for the reaction rate. Alternatively, the 
rate estimates could also be based directly on three-body 
resonance energies if available. 

An added benefit of the three-body formalism is that 
no assumptions regarding the preferred reaction path are 
needed. On the contrary, the most likely reaction mecha¬ 
nism can be deduced based on the structure of the angular 
wave functions corresponding to the lowest A„ eigenvalues. 
For most of the relevant energies the mechanism can be 
considered sequential through narrow two-body resonances. 
This is intuitively understandable if the three-body excited 
resonance state has a higher energy than the resonance(s) 
of the subsystem. On the other hand, a very low temper¬ 
ature only allows direct decay (or capture) since then the 
available energy is too low to populate the two-body res¬ 
onances even virtually. However, these reaction questions 
would eventually have to be answered by studying the dy¬ 
namic evolution from small to the large distances in more 
detail. 

In conclusion we have performed a full three-body anal¬ 
ysis of the nuclear structure and decay of critical waiting 
points in the rp process. This allowed us to study in detail 
the two-proton capture rate needed to bridge the waiting 


points, as well as the energy levels central in determining 
the effective lifetimes of these waiting points in a stellar 
environment. We find that a simple expression, based on 
an assumed Breit-Wigner shape, can accurately reproduce 
the full three-body rate calculation in the temperature re¬ 
gion 0.5 — 5 GK generally considered to be of astrophysi- 
cal interest. This led to a general method for estimating 
two-proton capture rates in the region around the critical 
waiting points based only on either the two- or three-body 
energies, if these are more readily available. Specifically, 
we predict, given the currently available experimental data, 
that the two-proton capture rate forming 70 Kr for temper¬ 
atures between 2 and 4 GK increase by a factor of 3.5 from 
0.4 • 10 _11 A^cm 6 s _1 mol -2 . 
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